An analytic model of the stereodynamics of rotationally inelastic molecular collisions 
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We develop an analytic model of vector correlations in rotationally inelastic atom-diatom collisions 
and test it against the much examined Ar-NO (X^II) system. Based on the Fraunhofer scattering of 
matter waves, the model furnishes complex scattering amplitudes needed to evaluate the polariza- 
tion moments characterizing the quantum stereodynamics. The analytic polarization moments are 
found to be in an excellent agreement with experimental results and with close-coupling calculations 
available at thermal energies. The model reveals that the stereodynamics is governed by diffraction 
from the repulsive core of the Ar-NO potential, which can be characterized by a single Legendre 
moment. 



Observing correlations among the vectors that charac- 
terize a collision can disclose all there is to know about 
how the collision proceeds [1]. Dudley Herschbach [2] 
likened vector correlations to "forbidden fruit" whose 
"tasting" reveals what would otherwise remain hidden. 
An example he frequently cites is the undoing of the az- 
imuthal averaging about the initial relative velocity vec- 
tor via a three-vector correlation, which reveals stereo- 
dynamical features lost by averaging over the initial dis- 
tribution of impact parameters. The pioneering work 
of Herschbach and coworkers [3] on vector correlations 
in the domain of molecular collisions spurred an effort 
to extract the information hidden in molecular dynam- 
ics computations, both quasiclassical and quantum, as 
these contain vector correlations as a default bonus [4]. 
However, even when characterized to the full by vec- 
tor correlations, the why of collision stereodynamics can 
only be answered as well as the theoretical method ap- 
plied to treat the collisions allows. In the present work, 
we extract vector correlations from an analytic model of 
direct rotationally inelastic atom-diatom collisions, and 
thereby gain a particularly simple, yet perspicacious in- 
sight into their stereodynamics. 

The collision model employed is based on the Fraun- 
hofer scattering of matter waves [5], recently extended 
to treat collisions in fields [6]. In contrast to classi- 
cal or semiclassical theories, the Fraunhofer model fur- 
nishes complex scattering amplitudes needed to extract 
the characteristics of vector correlations that reflect the 
quantum stereodynamics. Owing to its quantum nature, 
the model accounts for diffraction, interference, and other 
nonclassical effects. As an example, we deal with inelastic 
collisions of closed-shell atoms with rotationally polarized 
symmetric-top-equivalent linear molecules, represented 
by the much examined Ar-NO (X^U) systcm[7, 8]. The 
vector correlations obtained from the Fraunhofer model 
are found to be in an excellent agreement with the re- 
sults of experiments and close coupling calculations of 
Wade et al. [9]. This allows interpreting the collision 
stereodynamics of the Ar-NO (AT^H) system in terms of 
the Fraunhofer model. 

The Fraunhofer model was described in detail in 



refs. [5, 6]. It is based on the sudden approximation, 
which treats the rotational motion as frozen during the 
collision and thereby allows expressing the inelastic scat- 
tering amplitude in terms of the elastic one. The elastic 
scattering amplitude is, in turn, expressed in terms of 
the amplitude for Fraunhofer diffraction of matter waves 
from a sharp-edged, impenetrable obstacle acting in place 
of the molecular scatterer. At collision energies of hun- 
dreds of cm~^, consistent with the sudden approxima- 
tion, the shape of the scatterer is approximated by the 
repulsive core of the atom-molecule potential, with the 
attractive part disregarded. The Fraunhofer model ren- 
ders fully state- and energy-resolved scattering ampli- 
tudes and all the quantities that unfold from them in 
analytic form. 

The stereodynamics of an atom-diatom collision is 
usually described by a set of four vectors: the initial and 
final relative velocities, k and k', and the initial and final 
rotational angular momenta of the diatomic molecule, j 
and j'. We use the initial and final relative velocities k 
and k' to define the XZ plane of the space-fixed coordi- 
nate system, with the initial relative velocity k pointing 
along the Z axis. In keeping with the convention of Orr- 
Ewing and Zare [10], we characterize the spatial distribu- 
tion of the angular momenta relative to the XZ plane by 
the polarization moments ^^^^ , which arise as coefficients 
in the expansion of the density operator over the state 
multipoles [11]. Since within the Fraunhofer model the 
scatterer is two-dimensional, the model can only account 
for alignment, but not for orientation [6]. As a result, all 
polarization moments with odd fc or g vanish within the 
model. 

For this case study, we chose the k — k' — j' three- 
vector correlation in the Ar-NO {j ~ il = 1/2 ^ j' , O' = 
1/2) collisions, as this can be compared with the re- 
sults of experiments and close-coupling calculations of 
Wade et al. [9]. In addition, we illustrate the scope of 
the model by treating the k — j — k' and k — j — k' — j' 
correlations in the Ar-NO {j = n = 3/2 j', fi' = 3/2) 
scattering which, to date, have not been measured or 
evaluated. We restrict our considerations to the two low- 
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FIG. 1: Polarization moment Aq^'Q') pertaining to the 
k — k'— j' three-vector correlation in Ar-NO (j = fl = 
1/2, j', Q' = 1/2) collisions at 520±70 cm"\ The analytic 
results furnished by the Fraunhofer model (red solid line) are 
compared with the experiment (black dots) and close-coupling 
calculations (blue dashed line) of ref. [9]. 



FIG. 2: Polarization moment A2^(j') pertatining the 
k — k' — j' three-vector correlation in the Ar-NO {j — Q — 
1/2,^ = 1/2) collision at 520±70 cm"^ The analytic 

results furnished by the Fraunhofer model (red solid line) are 
compared with the experiment (black dots) and close-coupling 
calculations (blue dashed line) of ref. [9]. 



in ref. [9], and average over the e/f parity states as 
these have not been resolved in the experiment. We take 
into account the energy spread of the molecular beams, 
-E'coii = 520 ±70 cm~^, by averaging our results over three 
collision energies corresponding to the most probable en- 
ergy and to energies at half-maximum of an essentially 
Gaussian collision energy distribution. In determining 
the Ar-NO potential, we rely on the most recent po- 
tential energy surface (PES) obtained by Sumiyoshi et 
al. [12] and make use of only the average potential, Kjum, 
since the differential and depolarization cross sections are 
found to be only weakly affected by the difference PES, 
Vdif [13]. The PES of ref. [12] comes close to that of 
Alexander [8] and both PES's yield essentially the same 
polarization moments. 

TABLE I: Meaning of the A'~^^ (j) and A^l (j) aUgnment polar- 
ization moments. The Z axis points along the initial relative 
velocity k. The final relative velocity k' lies in the X > Q 
half of the XZ plane. The indicated ranges of the moments 
correspond to the high-j limit. The ^j-O) moment vanishes 
identically. 



Moment 






Meaning 
Range 


j along Z 
i±Z ^ -1 

i\\z ^ 2 


j along X or Y 

i\\x ^ -1 

i\\y ^ ^ 



In order to characterize the k — k' — j' three-vector 
correlation, we make use of the alignment moments 
Aq ' {y) and A2_|l(j') of the diatomic's final rotational an- 
gular momentum j' with respect to the XY plane. The 
Aq '{y) moment accounts for alignment of j' with respect 
to the initial relative velocity k and, in the high-j' limit, 
ranges between —1 and 2. Positive (negative) values of 

(2) 

Aq (j') correspond to j' ]| k (j' ± k, in which case j' lies 

(2) 

in the XY plane). The ^2+(j') nioment varies from — 1 
to 1 (in the high-j' limit). Its positive (negative) values 
correspond to alignment of j' along the F-axis (X-axis). 
This is summarized in Table I. 

(2) f2) 

Figures 1 and 2 display the Aq (j') and ^2+0') ™0" 
mcnts obtained in analytic form from the Fraunhofer 
model along with the results of experiment and close- 
coupling calculations of Wade et al. [9] . The agreement 
between the Fraunhofer model and the close-coupling cal- 
culation is compelling. 

Unfortunately, only two experimental points are avail- 
able for the scattering angles concerned. Therefore, it is 
not clear whether the oscillatory behavior at small-angles 
would indeed show up in an experiment. We hope that 
the present work will inspire an experiment whose reso- 
lution will suffice to clarify this issue. 

One can see that for zero scattering angle, i? = 0, 

(2) (2) 

Aq = —1 and A2_(^ = 0. The reason is geometric: in 
pure forward scattering, the j' vector must be perpendic- 
ular to k. Also, since k is roughly parallel to k' for very 
small "i?, the j' vector is approximately perpendicular to 
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FIG. 3: Polarization moment Aq^-'Q) pertaining to the 
k — j — k' three-vector correlation in Ar-NO {j = Q = 3/2, 

— 3/2) collisions at 520±70 cm~^ obtained from the 
Fraunhofer model. 




FIG. 4: The k — j — k' — j' four-vector correlation in the Ar- 
NO (j = n = 3/2, j', Q.' = 3/2) collisions at 520±70 cm~^ 
in terms of the dependence of the final polarization moment 
^o^'(j') on the initial alignment Aq'^'Q) = (black line), 
^o^'(j) = -0.8 (blue dashed line), and A[|^'(j) = 0.8 (red dot- 
ted line). Obtained analytically from the Fraunhofer model. 

k'. For small but nonzero scattering angles, i9 ~ 5°, the 

(2) 

Aq moment becomes positive, both for j' = 9/2 and 
17/2, indicating that j' tends to align along k. 

(2) 

The moment, on the other hand, exhibits narrow 
positive oscillations at very small angles (i? « 1°), but is 
in general negative, which corresponds to alignment of j' 



along the X-axis. 

Interestingly, the polarization moments presented in 
this paper are in a quantitative agreement with the ac- 
curate, close-coupling calculations, while the differen- 
tial cross sections for the Ar-NO scattering, evaluated 
in ref. [6], agree only qualitatively. From this we draw 
the conclusion that the polarization moments are mainly 
due to the hard-core part of the potential. This con- 
clusion is also supported by purely classical arguments 
based on the conservation of the projection of angular 
momentum on the collision's kinematic apse, see, e.g., 
ref. [14]. The phase shift of the moment's oscillations as 
derived from the model with respect to those obtained 
from the close-coupling calculation is likely due to ne- 
glecting in the model the attractive part of the Ar-NO 
potential. This explanation is supported by a general- 
ized Fraunhofer model that accounts for both attraction 
and repulsion [15] and which shifts the oscillations back 
toward smaller d. 

The alignment moments of Figs. 1 and 2 were normal- 
ized by the differential cross sections obtained from the 
Fraunhofer model. Since the Fraunhofer differential cross 
sections decrease faster with the scattering angle than 
their close-coupling counterparts [6], the oscillations of 
the analytic polarization moments are left relatively un- 
damped at large scattering angles. 

Since the j vector with j = 1/2 can only be oriented 
but not aligned, the only-to-alignment-sensitive Fraun- 
hofer model cannot handle vector correlations involving 
the j- vector with j = 1/2. However, since the j vec- 
tor with j ~ 3/2 can be aligned, we worked out the 
k — j — k' and k — j — k' — j' vector correlations for the 
Ar - NO {j = 3/2,0,^ j',n) scattering within the 

(2) 

il = 3/2 manifold. Figure 3 displays the (j) polariza- 
tion moment for the k — j — k' three- vector correlation. 
One can see that the small-angle scattering is favored by 

(2) 

positive values of the (j) moment, which corresponds 
to j II k, i.e., to a "broadside" approach of NO with re- 
spect to k, which enhances the scattering cross section. 

(2) 

However, for larger scattering angles, the Aq (j) moment 
becomes slightly negative, attesting to a preference for an 
"edge-on" approach with j -L k. 
Figure 4 exemplifies the k — j — k' — j' four- vector cor- 

(2) 

relation in terms of the alignment moment A^ '{j') of the 

(2) 

final j' for different polarizations (j) of the initial 
j. For an unpolarized initial state (black line), j' tends 
to align perpendicular to k ("broadside" recoil) for very 
small 7?, but reverses to a slight alignment in the parallel 
direction ("edge-on" recoil) for larger scattering angles. 
Initial polarization of NO such that j _L k is seen to result 
in only small changes of the final alignment (blue dashed 
line). However, in the case of a "broadside" approach, 
j II k, the stereodynamics changes significantly (red dot- 
ted line). The Aq (j') moment remains slightly negative 
throughout the range of scattering angles, indicating a 
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propensity for an "edge-on" recoil, with j' _L k. 

The Fraunhofer model readily explains the above re- 
sults: the analytic scattering amplitudes are proportional 
to the Bessel functions, which is a signature feature of 
diffraction. It is thus a diffractive oscillatory pattern that 
determines the angular dependence of the polarization 
moments. While the shape and frequency of the angular 
oscillations are entirely determined by the hard core of 
the PES, their position is somewhat influenced by the 
PES's attractive branch [15]. The Clebsch-Gordan co- 
efficients that appear in the scattering amplitude bring 
about selection rules that constrain the final parity of 
the states and the projections of the angular momentum 
j' on k. Within the model, the shape of the scatterer 
enters through the Legendre moments of a series expan- 
sion of the hard-core PES in terms of Legendre polyno- 
mials, Pk{cos6). The angular momentum algebra that 
the model entails gives rise to additional selection rules 
which allow for nonzero contributions to the polarization 
moments to arise only from Legendre moments of order 
K > j' — j- Therefore, the vector correlations for the 
j ~ 1/2, 3/2 j' ~ 9/2 channels are governed by the 
Legendre moment with k = 4, whereas the k = 8 Leg- 
endre moment governs the polarization moments of the 
j = 1/2, 3/2 j' = 17/2 channels. 

Moreover, since the Fraunhofer model can account for 
collisions in an electrostatic field [6], we investigated the 
effect of the field on the polarization moments. A field 
of 16 kV/cm, sufficient to significantly orient the NO 
molecule in the space-fixed frame, was found to cause 
only a tiny difference in the parity-resolved polarization 
moments as compared with the field-free ones. Upon 
averaging over the e// states, the effect of the field was 
found to be altogether negligible. 

In summary, we made use of the Fraunhofer model 
of direct rotationally inelastic atom-diatom collisions to 
study vector correlations in such collisions analytically. 
The vector correlations obtained from the model closely 
reproduce those extracted from close-coupling calcula- 
tions which, in turn, agree well with experiment. The 
Fraunhofer model of vector correlations demonstrates 
that the stereodynamics of the Ar-NO rotationally in- 
elastic collisions is contained solely in the diffractive part 
of the scattering amplitude which is governed by a single 
Legendre moment characterizing the anisotropy of the 
hard-core part of the system's PES. Given the "geomet- 
ric" origin of this behavior - ordained by the angular 
momentum algebra ~ we expect to find a similar behav- 
ior in other systems. 
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